On the determination of curvature and dynamical Dark Energy 



J.-M. Virey 1 , D. Talon- Esmieu 2 , A. Ealet 2 , P. Taxil 1 and A. Tilquin 2 

1 Centre de Physique Theorique*, Case 907, F-13288 Marseille Cedex 9, France 
and Universite de Provence, Marseille, France 

2 Centre de Physique des Particules de Marseille, Case 902, 
F-13288 Marseille Cedex 9, France 



Abstract 

Constraining simultaneously the Dark Energy(DE) equation of state and the curvature 
of the Universe is difficult due to strong degeneracies. To circumvent this problem when 
analyzing data it is usual to assume flatness to constrain DE, or conversely, to assume 
that DE is a cosmological constant to constrain curvature. In this paper, we quantify the 
impact of such assumptions in view of future large surveys. We simulate future data for 
type la Supernovae (SNIa), Cosmic Microwave Background (CMB) and Baryon Acoustic 
Oscillations (BAO) for a large range of fiducial cosmologies allowing a small spatial cur- 
vature. We take into account a possible time evolution of DE through a parameterized 
equation of state : w(a) = wo + (1 — a)w a . We then fit the simulated data with a wrong 
assumption on the curvature or on the DE parameters. 

For a fiducial ACDM cosmology, if flatness is incorrectly assumed in the fit and if the 
true curvature is within the ranges 0.01 < Qk < 0.03 and —0.07 < Qk < —0.01, one will 
conclude erroneously to the presence of an evolving DE, even with high statistics. 
On the other hand, models with curvature and dynamical DE can be confused with a flat 
ACDM model when the fit ignores a possible DE evolution. We find that, in the future, 
with high statistics, such risks of confusion should be limited, but they are still possible, 
and biases on the cosmological parameters might be important. 

We conclude on recalling that, in the future, it will be mandatory to perform some com- 
plete multi-probes analyses, leaving the DE parameters as well as the curvature as free 
parameters. 
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1 Introduction 

Since the discovery of the acceleration of the Universe [1, 2] many theoretical interpreta- 
tions have been developed to find an explanation. These models introduced in general a 
new component called Dark Energy (DE) with an unknown nature [3, 4]. The distinction 
between the various interpretations is critical for future cosmology and for its connections 
with fundamental physics. 

Many studies concentrate on the equation of state (EoS) w, defined as the pressure 
to density ratio of this new component. This parameter is equal to —1 for a classical 
cosmological constant but can be different and/or can vary with cosmic time depending 
on the various DE models (for a review see e.g. [4]). On the experimental side, the 
determination of the nature of the Dark Energy by constraining w(z) is largely debated 
and should use a large variety of probes and methods [5]. 

The WMAP data combined with the recent SDSS data [6, 7, 8, 9] are currently in 
very good agreement with an adiabatic ACDM model. However, most of studies in the 
literature assume a spatial flatness to constrain the dynamical Dark Energy parameters 
or, conversely, assume a cosmological constant to derive constraints on curvature. At 
the current level of statistical precision, this can be justified but various authors have 
emphasized the dangers one can encounter in these too simple strategies [10, 11, 12] and 
the impact of possible biases, when the statistical errors start to be small, should be 
investigated. 

In general, extracting simultaneously w(z) and the curvature (and also the matter 
density) is not possible since these parameters are perfectly degenerated. This is often 
referred to the so-called "geometrical degeneracy" [13, 14, 15, 16]. As almost all in- 
flationary scenarios predict spatial flatness < 10~ 10 ), this property is often taken 
for granted in the literature. However the precision is not so compelling. Experimen- 
tally, the error on the curvature parameter as given by the WMAP CMB data, 
is ±0.006(68%); ±0.013(95%) when adding a Supernovae sample, the Baryon Acoustic 
Oscillations (BAO) from SDSS and assuming a cosmological constant [9]. Assuming a 
parameterized w(a) = w + (1 — a)w a equation of state (see below), the precision is only 
of ±0.018(68%); ±0.04(95%) [9, 17]. Relaxing on the combination with SDSS and the 
Supernovae's (SNIa), leads even to a precision of ±0.10 [9]. Then, the precision on cosmic 
curvature is very dependent on the assumption concerning Dark Energy, and is strongly 
limited today by the number of possible free parameters. 

On the other hand, concerning Dark Energy evaluation, even the more recent standard 
analysis as the WMAP ones [8, 9] or SNLS [18], give results on the Dark Energy parameters 
using the flatness hypothesis and/or no Dark Energy evolution. These analysis use, as is 
now common, the compressed parameters : A from BAO [6] and/or of the shift parameter 
R to include the WMAP CMB data [8, 9, 17, 19]. 

Few recent studies [20, 21, 22] have extracted simultaneously a free curvature and a 
dynamical Dark Energy from the present data. Such combined analyses need always some 
external priors or/and assumptions (prior on h, the present value of the Hubble constant in 
unit 100 km. s~ 1 .Mpc~ 1 , assuming massless neutrinos, adiabatic perturbations only...). It 
is known for some time [23] that using incorrect priors can lead to an artificial convergence 
toward a particular model, e.g. the ACDM model, (for the consequence of using a strong 
prior on fl M , see e.g. [24, 25] and [26] for a recent contribution). The situation could be 
worse when the number of relevant parameters is increased e.g. when relaxing a constraint 
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on the curvature and/or the nature of Dark Energy. 

In this paper we adopt the same kind of attitude as in [24, 27], focusing on the question 
of the interplay between the curvature and some DE parameters. We extend the work of 
previous authors (e.g. [11, 12]), by addressing these questions more quantitatively and 
taking into account more realistic observational conditions. 

To do parameter bias estimation with free curvature, which means to compare the 
fitted parameters values to fiducial ones, a combined analysis is mandatory. We choose to 
stay in a frame of three combined probes : type la Supernovae (SNIa), Cosmic Microwave 
Background data (CMB) and Baryon Acoustic Oscillations data (BAO). This combina- 
tion (which is independent of H ) has enough statistical power to start to relax on the 
flatness hypothesis. Furthermore, we take advantage on the compressed parameters used 
in the literature as they contain most of the cosmological information needed to reduce 
Dark Energy parameters errors at the same level as from a complete analysis [17]. In a 
simulation study, where all other parameters are fixed, this will not introduce additional 
biases with the main advantage of avoiding a complete CMB and BAO data analysis. 
Therefore, we use the R shift parameter [13, 15, 28, 17] for the CMB, and the so-called 
A parameter for the BAO [6]. This will reduce drastically the computing time, which 
is useful since we want to simulate a large range of models, and it will not affect the 
generality of our analysis. 

In section 2, we present the simulated sample of SNIa, the R and A parameters, and 
the statistical method used for minimization. We also introduce a detectability criterion, 
based on the goodness of the \ 2 value of the fit. 

We study in section 3 and 4 two cases which illustrate the limit of some existing 
analyses. First, we assume flatness to constrain a dynamical Dark Energy in a case of an 
Universe with a non flat curvature and a cosmological constant. Second, we assume that 
Dark Energy is a cosmological constant to constrain curvature in a model of dynamical 
Dark Energy. This case is even more complicated since the number of fiducial parameters 
has increased compared to the previous scenario, this is the "triple trouble" situation 
mentioned in [11]. 

We conclude in section 5 that to get any strong conclusion on Dark Energy evolution 
one will need to fit simultaneously curvature and dynamical Dark Energy parameters in 
a multi-probe analyses. 

2 Cosmological probes, fitting method 
2.1 The Dark Energy parameterization 

To simulate and fit the involving DE models, we need a parameterization of the EoS 
w(z). We should choose for our study a parameterization that takes into account a 
possible evolution on all redshift range up to the CMB, with a minimum number of free 
parameters to avoid too large errors. Indeed it has been shown that at most a two- 
parameter model can optimally be constrained by future data [29, 25]. We choose the 
Chevallier-Polarski-Linder (CPL) parameterization [30, 31], which has been extensively 
used in the literature, in particular by the Dark Energy Task Force (DETF) [5], as a 
phenomenological benchmark to compare and contrast the performances of different DE 
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probes (see e.g. [32]): 

w(z) =w + w a z/ (l + z)=w + (l- a)w a (1) 

a being the scale factor. The free parameters w and w a characterize the DE model. The 
cosmo logical constant corresponds to w = — 1 and w a = 0. 

Despite its simplicity, this parameterization, which is bounded for high redshifts 
(w(z — > oo) = wq + w a ), fits well the evolution predicted by a large class of quintessence 
models (see [33] and also [34] for a recent discussion). Note that the CMB analysis 
introduces in general a high redshift constraint for the chosen Dark Energy parameteriza- 
tion [25], namely w(z C mb) < at last scattering. With our parameterization we should 
impose Wq + w a < (we refer in the following to this limit as the " CMB boundary con- 
dition"). Note however that some DE models can be characterized by an EoS such that 
w(zcmb) > 0. This is the case of the early Dark Energy models characterized by a non- 
negligeable DE density at last scattering [35] and also for the tracker quintessence models 
[36] which have null and positive EoS at early time as shown explicitly in [37]. Therefore 
we will allow to consider sometimes the part of the parameter space where the "CMB 
boundary condition" is violated. 

Concerning our simulations, this particular choice of the parameterization has no 
strong impact on the result as soon as fiducial models and fitted parameters are using the 
same modelisation. The use of another two-parameter parameterization will not change 
the main message of this work, we have checked that it can only slightly modify the size 
of the errors. 



2.2 SNIa samples 

In the standard Friedmann-Roberston- Walker metric, the apparent magnitude of astro- 
physical objects can be expressed as a function of the luminosity distance: 



m(z) 



5 log 10 (D L ) + M B - 5 log w (H /c) + 25 = M s + 5 log 10 (D L ) 



(2) 



where Mb is the absolute magnitude of SNIa, M s may be considered as a normalization 
parameter and Dl(z) = (H /c) cIl(z) is the H -independent luminosity distance to an 
object at redshift z. It is given by: 



VP. 

where S(x) = sinh(x), x, sin(x) for Qk > 0, =0, < respectively, with 

= 1 — VLm — &x 

Qm and fix being respectively the reduced matter and DE densities. We have: 



E(zf 



where 



H(z] 
Ho 

Px{z) 
Px(0) 



'i + z yn M + 



= exp 



Px(0) 

3 / (l + w(z')) d ln(l + 
. Jo 



n x + (i + zyn k , 



(3) 



(4) 



(5) 



(6) 



4 



Note that we have neglected the radiation component and this has no impact on 
our results. With the parameterization given in eq.(l) the last equation becomes: 

Px( Z ) — ( i _i_ _\ 3(l+wo +w a ) ~3w a z/(l+z) ( 7 \ 

Px{0) " { ] ' { ) 

The studies presented in this paper are performed with simulated Supernovae samples 
with statistics equivalent to what we expect to have in the future. We concentrate on two 
sets of data which simulate the statistical power of the forthcoming and future data. 

• We simulate forthcoming data from ground survey as the large SNLS survey at 
CFHT [38, 18]. This survey has started in 2003 and the estimation after 5 years 
of running is to register a sample of 700 identified SNIa in the redshift range 0.3 < 
z < 1. We simulate a sample as described in [27]. This corresponds to our "short 
term" scenario. 

• We simulate data from a future space mission like JDEM/SNAP [39], which plans 
to discover around 2000 identified SNIa, at redshift 0.2< z <1.7 with very precise 
photometry and spectroscopy. The Supernovae distribution is given in [40]. This 
corresponds to our "long term" scenario. 

The magnitude dispersion is assumed to be constant and independent of the redshift 
at the level of 0.15 for all Supernovae after correction, for both samples. We have also 
neglected systematical errors in this study. 

A set of very well calibrated SNIa at redshift < 0.1 should be measured by the SN 
factory collaboration [41]. This sample is needed to normalize the Hubble diagram and 
will be called in the following the " Nearby" sample. A " SNAP" (" SNLS" ) sample means 
in this paper a simulation of the statistics expected from the JDEM/SNAP like mission 
("SNLS" survey) combined with the 300 (200) nearby SNIa expected to be measured at 
the time of SNAP (SNLS). 



2.3 The CMB and BAO constraints 

We need to add simulated CMB data to take advantage of the current and future con- 
straints on curvature. We have simulated the CMB constraints by using the CMB shift 
parameter R [13, 15, 28, 17] which is the scaled distance to recombination: 

where zqmb = 1089 is the redshift at the epoch of recombination. This parameter contains 
the full geometrical degeneracy inherent to the CMB [15]. It is measured to be Rwmap = 
1.71 ± 0.019 from the 5-years WMAP data [9]. 

Using R we lose some pieces of information which are potentially present in the CMB 
but it contains enough statistical power to constrain the Dark Energy sector at a level 
comparable to the one of a full CMB analysis. To evaluate the relevance of our results 
and of this assumption, we have made some full CMB calculations and compared the 
results when evaluating the Dark Energy parameters in some models. We found a correct 
agreement between the two methods. For a comparison between the use of the CMB R 
parameter and full CMB calculations, one can refer to [28, 17]. 
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We consider BAO measurements to add an additional constraint on the matter den- 
sity Q M . We use the so-called A parameter, described in [6], which encodes almost all 
information on the cosmological parameters from present BAO observations : 



A = 



E(z BAO y/ 3 



i 



ZBAO 



2/3 



(9) 



where z BAO = 0.35. It is measured as A SDS s = 0.469(n s /0.98)~ a35 ±0.017 with n s = 0.95 
[6]. The definition of A is an approximation which incorporates a mixture of the trans- 
verse and radial information present in the data. 



In practice we calculate the values of R and A for each simulated fiducial model. 
These values are used in the fits along with the errors expected for each observational 
scenario. For the mid term scenario we take the present errors on R and A: u{R) = 
0.019 (WMAP)[9] and a(A) = 0.017 (SDSS) [6]. For the long term scenario we choose 
<t(R) = 0.01 as roughly expected [42] from the Planck mission [43], and cr(A) = 0.005 as 
expected for future large surveys like a " stage- IV" space based mission dedicated to BAO 
as described by the DETF [5, 44]. 



2.4 The statistical method 

To proceed in practice, we choose a fiducial model that is a set of fiducial parameters 
and simulate data (SNIa magnitudes, R and A). In the most general case the fiducial 
parameters are Q^-, Q^, Wq and w%. Note that the value of the SNIa normalisation 
parameter Mf is taken to -3.6: its precise value is not much relevant for the following. 
For definiteness we will assume Vt^ = 0.3 throughout the paper except otherwise specified 
but we have checked that reasonable deviations (±0.2) from this value has a very small 
influence on our conclusions. 

Then, having chosen a fiducial model, we simulate data for this model and we fit 
the simulated data with a wrong assumption on the curvature or the DE parameters. 
Finally, we compare the fiducial and fitted values of the various parameters along with 
the respective errors, to estimate how much our conclusions could be affected by the wrong 
assumption. 

To analyse the (simulated) data, a minimization procedure has been used [45]. A 
standard Fisher matrix approach allows a fast estimate of the parameter errors. This 
method is however limited as it does not yield the central values of the fitted param- 
eters. Then we adopt, unless specified, a minimization procedure based on a least 
square method. The least square estimators are determined by the minimum of the 
X 2 = (m — M(z, f2, u>) T V _1 (m — M(z,Q,w)), where m = (m 1 ...m n ) is a vector which 
contains the simulated magnitudes plus R and A. The vector M(z, Q, w) is the corre- 
sponding vector of values for the fitted parameters and V the covariance matrix. The 
errors on the cosmological parameters are estimated at the minimum by using the first 
order error propagation technique: U = J.V.J T where U is the error matrix on the 
cosmological parameters and J the Jacobian of the transformation. A full n-parameter 
fit (n-fit) of the simulated data gives the central values and errors for the n parameters, 
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along with their correlations which are in general strong. 



2.5 Detectability criterion 

We want to identify among models the more problematic ones where the fit fails to identify 
a problem. Indeed, in any data analysis, a wrong assumption can be early detected 
through a simple x 2 test: a high x 2 indicates that the goodness of the fit is bad. In our 
context, when fitting simulated data with some hypotheses, a high x 2 wm indicate that 
some of these hypotheses are wrong. 

However, rejecting models from a x 2 test with real data is always tricky, therefore we 
demand a high level of detectability. Consequently, models where the x 2 test gives a 5<r 
effect, (x 2 > 5o"(x 2 )) where the rms of the x 2 is °"(x 2 ) — a/2^ / and where N^of is the 
number of degrees of freedom in the fit, are considered to be easily identified as in conflict 
with the input hypothese. 

On the other hand some models passing this test can be biased as the fitted parameters 
are too far from the fiducial ones. In the next sections, we identify the models which yield 
a good value of the x 2 even if the fitting hypotheses are incorrect. A careful attention is 
given to the resulting confusions and biases. 

3 Reconstruction of DE dynamics with a wrong as- 
sumption on curvature 

If we consider realistic experimental conditions, the reconstruction of an evolving Dark 
Energy is difficult when we leave the curvature free due to a complete degeneracy between 
the relevant parameters. This comes from the integral relation between cosmological 
distances and the cosmological parameters [23] and the only way to solve this question 
is to fix some parameters or/and to combine different informations. However, such a 
strategy can lead to some biased interpretations. This has recently been studied on 
theoretical grounds in [12]. Following the same approach, we simulate fiducial models as 
ACDM Universes with non-zero curvature and we calculate the errors on the cosmological 
parameters as expected in future observational conditions. In this way we are able to 
quantify when some models with curvature, consistent with the CMB constraint, can 
start to bias the results on the DE EoS when flatness is incorrectly assumed in the fit. In 
practice, we fix in the range 0.1 — 0.5 (illustrations are given for Qfj = 0.3) and we 
vary to scan f2f between -0.1 and 0.1. Then we perform the fit of Qm,u! and w a , 
assuming flatness. 

To test the wrong hypothesis we first test the goodness of the fit as described in the 
previous section. When fiducial curvature is too far from flatness, the x 2 is too high 
and the wrong assumption is detectable. Models with Q k lower than -0.08 (-0.07) and f2 fc 
greater than 0.05 (0.03) for the short term (long term) scenario are detected by the simple 
X 2 test. We study in more details some models which give good x 2 fits to determine the 
risks of a wrong conclusion on the values of the DE parameters. 

In Figure 1 we present 2 examples of fits with good x 2 m the (wq, w a ) plane. We have 
chosen two different fiducial models: a ACDM open model (upper plots) with (wq — — 1, 
w f = 0, Vt F M = 0.3, tt F = 0.68) and a ACDM closed model (bottom plots) with (wg = 
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Figure 1: Illustration of the bias problem in the (wq, w a ) plane for a slightly open model 
(upper plots) with Qf^ = 0.3, Q x = 0.68 = 0.98) and for a slightly closed model 
(lower plots) with Qf^ = 0.3, Q x = 0.72 (Q^ = 1.02). Results are given for a short term 
scenario (left plots) and a long term scenario (right plots). A fiducial ACDM model has 
been assumed. Contours are at la and correspond to some fits with (solid) or without 
(dashed) the flatness assumption. 

— l,wj[ = 0, = 0.3, Q x = 0.72). In each case the full simulated data (SNIa, R and 
A), are fitted with the (wrong) hypothesis of flatness. The la contours correspond to the 
solid curves and the left (right) plots correspond to the short term (long term) scenario. 

One can see that the fitted w and w a central values are far from the fiducial values in 
each case. We get for positive curvature (Qk < 0) a wq value which is pushed towards the 
phantom regime wq < — 1 and a w a value which is positive, and the converse for negative 
curvature (f2 fc > 0) i.e. that w > — 1 and w a < 0. Let us stress again that this bias is 
not detectable with true data since the fit is good (low y 2 ). Note that in the same time 
the Qm value is well reconstructed. 

For the long term scenario, the statistical power increases, then the error ellipses are 
obviously smaller but the bias is higher. For comparison we show the contours resulting 
from a complete fit where the curvature is a free parameter (dashed contours). These 
contours are centered on the fiducial values as expected. In the long term scenario the 
contours from the two approaches are clearly disconnected, hence the problem becomes 
stronger. 

Let us quantify more precisely the biases which are expected on these two DE EoS 
parameters. We vary the curvature and determine the zones of bias, defined as | Wq + 1 | > 
a(w ) and | w a |> a{w a ). The DE parameters bias means that a cosmological constant 
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Bad Fit 


Biased Zone 


Validity Zone 


Biased Zone 


Bad fit 


short term 


ft£ <-0.08 


-0.08< Sl% <-0.02 


-0.02< Q% <0.02 


0.02< n b k <0.05 


Q% >0.05 


long term 


ft£ <-0.07 


-0.07< n% <-0.01 


-o.oi< n£ <0.01 


0.01< <0.03 


ft£ >0.03 



Table 1: Validity range of the flatness assumption for the determination of the DE EoS 
parameters wo and w a . 

with some amount of curvature will be interpreted as a dynamical DE model in a flat 
Universe. The results are summarized in Table 1: 

• In the positive curvature sector, —0.08 < Q k < —0.02 (short term), —0.07 < Q k < 
—0.01 (long term). The upper bound is the bias limit, above this value the results 
are not biased. The lower bound is the goodness of the fit limit, below this value 
the x 2 explodes. 

• In the negative curvature sector, 0.02 < Q k < 0.05 (short term), 0.01 < f2f < 0.03 
(long term): now the lower bound is the bias limit and the upper bound the goodness 
of the fit limit. 

In other words, if we do not take into account the current and future curvature uncer- 
tainties, results on the DE parameters can be biased even in a long term scenario. This 
implies a potential misinterpretation of the DE dynamics even when using a combination 
of precise probes both at low and high redshift, as soon as a spatial flatness is assumed. 

We conclude that w(z) can be reliably measured with the flatness assumption only if 
the true curvature is \Q k \ < 0.02 (< 0.01) for a short (long) term scenario. 

4 Reconstruction of curvature with a wrong assump- 
tion on DE 

Another common strategy is to pin down the curvature of the Universe when assuming 
that DE is a cosmological constant (ACDM hypothesis). In this section we will investigate 
the biases introduced by this assumption. 

After a short illustration of the problem, we study in details the possible confusion of 
some non-flat dynamical DE models with a flat cosmological constant. Then we quantify 
the models which can provide a biased estimation of the cosmological parameters in an 
undetectable way. 



4.1 Illustration 

In this study, we test the ACDM model assumption by a fit of Qm and fix- First we 
perform the fit for two examples which illustrate the situation. We consider two models 
with a dynamical DE for two different curvatures (f2f = 0.03 and Q k = —0.02), still 
fixing the matter density to = 0.3. 
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We show on Figure 2 the results of the fits in the (Qm, ^x) plane for the combination 
of SNIa, R et A in the short term scenario. The straight line corresponds to flat models 
(Q M + Q x = 1)- We show in solid curves the contours obtained from a fit assuming that 
DE is described by a ACDM model. We see that the fitted values of the matter and 
Dark Energy densities (crosses) are far from the fiducial ones (stars) due to the wrong 
hypothesis. Therefore the fitted values of curvature are compatible with the flatness at 
la. For comparison the fits without assumption on the DE are performed and the results 
are plotted with the dotted contours. Note that these ellipses are completly disconnected. 
These results are given for a short term scenario but this bias still exists for the long term 
scenario. 




Figure 2: Contours at la in the (fyw, &x) plane for a short term scenario. The left (right) 
plot corresponds to the fiducial cosmology Vt^j = 0.3, Q x = 0.67, that is f2f = 0.03, 
w f = _0.9 and iu£ = 0.5 = 0.3, Vt x = 0.72, that is = -0.02, w£ = -0.9 and 
w% = —1.5). Plain (dotted) curves are obtained from fits with (without) the assumption 
that DE is a cosmo logical constant. The straight line corresponds to flat universes. 

With this illustration we have shown that it is possible to interpret a dynamical DE 
model in a Universe with a non-zero spatial curvature as a flat cosmological constant 
model, for present statistics as well as for future surveys. In the following we look for 
the DE models and the possible amount of curvature where such a misinterpretation of 
combined data is possible. 



4.2 Confusion with flatness 

We want to identify the DE models with non zero curvature which can be compatible 
with a flat ACDM model, when assuming incorrectly w = — 1 in the fit. 

For a given set of fiducial parameters (fif, Wq, w£) we perform a fit imposing w — — 1. 
For each tested model we first verify the goodness of the fit and then we determine if the 
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fitted value of the curvature is compatible with the flatness by satisfying \Q k \ < o-(Q k ). 

We present in Figure 3 the results for four fiducial curvatures (Q k = 0.06, 0.02, —0.02, —0.03) 
in the short term scenario. Models with a good fit are located between the two solid curves 
and represent a large area of the plane. It means that many DE models in a curved space 
may give a correct fit for a cosmological constant. For some of these models it exists 
a confusion with a flat Universe, i.e. these non-flat models are interpreted as flat ones. 
For each fiducial curvature these models form a confusion zone in the plane. Note that 
the models chosen for the illustration of the problem on Figure 2 are marked by stars on 
Figure 3 and are within the confusion zones. 




Figure 3: The black contour gives the limit of goodness of the fits. Other contours 
correspond to the zones of confusion with a flat ACDM model for the following fiducial 
cosmologies = 0.3): fi£ = 0.64 (Qf = 0.06) (plain), fi£ = 0.68 (fif = 0.02)(dotted- 
dashed), fi£ = 0.72 (fif = -0.02) (long-dashed) and fi£ = 0.73 (fif = -0.03) (short- 
dashed). The straight line is the "CMB boundary condition" w + w a < 0. Calculations 
are done for the short term scenario. 

We have performed the same study for the long term scenario. The results for fiducial 
curvature (fif) between -0.01 and 0.03 are given in Figure 4. Thanks to the higher 
statistics, the area of the (wq,w^) plane where the fit passes the x 2 test is smaller. 
Therefore the wrong assumption is detectable for a greater number of models. 

Concerning models which pass the x 2 test but may be confused with a flat cosmological 
constant model, we remark that DE models with positive curvature (fl k < 0) are in the 
lower part of the plane (roughly Wq > — 1 and w% < 0) on both Figures 3 and 4. With 
high statistics, most of these problematic models exit the zone of goodness of the fit and 
can be rejected as soon as Q k < —0.01. 

The DE models with negative curvature (Q k > 0) that may be confused with a flat 
cosmological constant model are in the upper part of the plane (roughly Wq < — 1 and 
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Figure 4: The black contour gives the limit of goodness of the fits. Other contours 
correspond to the zones of confusion with a flat ACDM model for the following fiducial 
cosmologies {VL F M = 0.3): fi£ = 0.67 (fif = 0.03) (plain), fi£ = 0.69 = 0.01)(dotted- 
dashed) and fi£ = 0.71 = -0.01) (dashed). The straight line is the "CMB boundary 
condition" w + w a < 0. Calculations are done for the long term scenario. 

w% > 0). Note that some confusion zones are close to, or above, the "CMB boundary 
condition". If we ignore this limit, some very open models (up to Q F = 0.6) can still give 
a correct x 2 an d be confused with a flat cosmological constant Universe. However, for 
these models, the size of the confusion zone is very tiny and corresponds to a very small 
phase space volume. 

Therefore, our analysis cannot exclude that some open dynamical DE models with 
a large departure from flatness can induce a confusion with a flat ACDM model. The 
proximity or violation of the CMB boundary condition indicates that these DE models 
are rather exotic: of both phantom and early Dark Energy types. The smallness of the 
phase space volume for these models show also that these exotic DE models should be 
fine tuned to give a confusion. This situation means that our analysis has reached its 
limits and that a better description of the DE model (beyond the CPL parameterization) 
and of the CMB (beyond the use of R) should be done to yield a confirmation of these 
possible confusions. 

One can find in the literature constraints on curvature coming from different combined 
analyses where the DE dynamics is left free [20, 21, 22]. They give Qk < 0.02 at 95%CL, 
concluding that open models are disfavored. Our last result indicates that one has to be 
cautious: constraints depend of the statistical methods which are very sensitive to the 
size of the phase space volume [25] and the current combined analyses of data may have 
neglected or missed these open exotic DE models. 
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In summary, even with a large statistical sample, it is still possible to confuse a dy- 
namical DE non-flat model with the flat ACDM model when a cosmological constant is 
assumed. This is due to the geometrical degeneracies between the DE and the CMB 
parameters and this cannot be completely solved by the use of a combined analysis with 
SNIa and BAO, even at high statistics. To draw any conclusion on the description of 
Dark Energy, it is mandatory to fit enough DE parameters with the curvature as a free 
parameter. 

4.3 Bias on the reconstruction of the reduced densities 

In this section, we study the biases on the reduced densities and Q x , still in the context 
of fiducial curved dynamical DE models analyzed assuming a cosmological constant. For 
that purpose, we vary in the range 0.1 to 0.5 and Q x in a range such that | < 0.2. 
Then, we scan the (wq,w f ) plane. We quantify the biases introduced by the wrong 
hypothesis, when trying to reconstruct the cosmological parameters Qm and Qx- We 
define the bias of a parameter p by B p = \p F — p\ and we say that p is biased (valid) if 
the bias is larger (smaller) than the error obtained for p i.e. if B p > cr(p) (B p < u{p)). 
(One can consult [27] for more details on these definitions). 

There are some zones where the x 2 is too high, then the fit is not good and the problem 
is detectable. When the x 2 is good we can define two types of zone : 

- the zone where the x 2 is good and the fitted parameters are well reconstructed even 
with the assumption of w = — 1, we call it the Validity Zone (VZ). This zone is centered 
on the A point (w F = — 1, w F = 0) where the fit assumption is true. 

- the zone where the x 2 is good but all fitted parameters are wrong 1 . The previously 
discussed zones of confusion belong to this zone and correspond to the cases where we 
obtain erroneously Qk = 0. 

Figure 5 shows one example where = 0.25 and fl x = 0.7 in the long term scenario. 
The result is generic for all models considered in this study and whatever the statistic is. 
We see that, with the combined analysis and high statistics, there are large zones in the 
plane where the cosmological parameters are wrong, if we assume DE to be a cosmological 
constant. This fact is well known. The new information is that many models will be 
interpreted as a ACDM because of the degeneracy of curvature and dynamical Dark 
Energy and this will not be detectable. 

The position of the various zones in the fiducial plane can move when fiducial param- 
eters like f2£f or/and Q x are changing: 

- An increase of £l x reinforces the power of the x 2 test (i.e. the goodness of fit zone 
decreases) but reinforces also the bias problem (i.e. the VZ decreases). 

- Changing Q 1 ^ shifts the location of the various zones in the fiducial plane. 

Scanning all possibilities, we can always find some values for the densities where the 
X 2 is low, and the fitted results are valid or biased. Our study shows also that the addition 
of SNIa and BAO reduces the degeneracies but do not eliminate them. 

1 For the short term scenario another situation may appear where fife is well-reconstructed but Hm and 
fix are biased, namely the biases compensate to give accidentally the correct curvature. This possibility 
disappears at high statistics. 
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Figure 5: The black contour gives the limit of goodness of the fits. Other contours 
correspond to the zones of confusion with a flat ACDM model (labelled by CZ), the validity 
zone (labelled by VZ) and the biased zone (labelled by BZ). The fiducial cosmology is 
Vt F M = 0.25 and fi£ = 0.7 = 0.05). The straight line is the "CMB boundary condition" 
Wo + w a < 0. Calculations are done for the long term scenario. 

Finally, let us notice that we find also that assuming only a constant w, different from 
— 1, leads to similar confusion and biases. 

5 Conclusions 

The vicious circle of assuming flatness to constrain dark energy, or conversely, of assuming 
that dark energy is a cosmological constant (or have a constant EoS) to constrain cur- 
vature, should be avoided. Indeed, such assumptions may provide very important biases 
and confusion if they are not true. 

We have used a combined analysis using SNIa, BAO and CMB data, taking avantage 
of the use of the compressed parameter A for BAO and shift parameter R for CMB. The 
combination will help to break degeneracies, origin of the biases, which come mainly from 
the integral relation between cosmological distances and the cosmological parameters [23] . 

We see that assuming flatness, a cosmological constant model with a small amount of 
curvature can be interpreted as a dynamical DE model in a flat universe. We find that, 
even using a multi-probe approach and with future high statistics, such hypothesis will 
bias the interpretation of the data if the true curvature is in the range 0.01 < Q k < 0.03 
or —0.07 < Q k < —0.01. This is due to the fact that, even with small statistical errors, 
this range of non-flat models will give a good x 2 when assuming flatness, because of the 




14 



degeneracy with the DE sector. We emphasize that these models are always inside the 
current level of uncertainties of WMAP data, (when an evolving DE is assumed in the 
analysis). Let us also state that the PLANCK mission will not allow to reduce strongly 
these uncertainties as the information on the curvature is mainly coming from the first 
peak which is already quite precisely determined. 

Conversely, if we suppose that Dark Energy is a cosmological constant, a dynamical 
non flat DE model may be interpreted as a flat cosmological constant model. We show 
that for closed models, there is no confusion as soon as £l k > —0.01. For open models, 
there is a greater risk of confusion, independently of the statistical sample, in particular 
if 0.01 < Qk < 0.05. However, the confusion may still happen at > 0.05, but only for 
very exotic DE models (of both phantom and early DE types) which seem less favored by 
the current data. 

In addition, we show that in case of misinterpretation, when assuming a cosmological 
constant, the cosmological parameters Vt M and Q x will be biased and this will be un- 
detectable in the fit for a large range of models. In any case, high statistics (long term 
scenario) help, in the sense that they allow to detect more easily a wrong hypothesis 
through the x 2 test. However, there is still a number of DE models where biases will not 
be detected and the smaller errors would increase our confidence in a wrong interpretation. 

In general, in most of the analysis of real data, two fits are performed independently : 
one with the assumption of flatness and the other with the assumption of a Cosmological 
constant. Consistency between the two fits is then considered as a proof a reliability. 
However, this approach is not exempt of drawbacks. For example, in the spirit of our 
analyses of Section 3, if we vary the DE parameters beyond a simple cosmological con- 
stant we find some confusions zones which overlap the ones of Figures 3 and 4. This 
means that these models yield a good \ 2 m & U cases and point toward a flat ACDM. In 
these cases the circular logic [12, 10] of the "consistency" test fails. 

The strategy which consists of relaxing the curvature and the DE parameters within 
one global procedure will be the only one which avoids confusion. In addition, when 
comparing the errors expected on the cosmological parameters in the fit procedures with 
or without the wrong assumptions (see the plain and dashed contours of Fig. 1 and 
Fig. 2), one sees that in a combined analysis using SNIa+CMB+BAO, there is no gain of 
assuming flatness or/and that DE is a cosmological constant. This fact has been already 
noticed, see e.g. [11]. 

In conclusion, our quantitative analysis allows us to strongly support the points already 
stressed by a few authors : it is imperative to let the curvature and the DE parameters as 
free parameters. This is already true for the analysis of present data as well as for data 
expected from the future large surveys. 
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